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Abstract. We report the TeV gamma-ray observations of the nearby normal spiral galaxy NGC 253. At a distance 
of ~2.5 Mpc, NGC 253 is one of the nearest starburst galaxies. This relative closeness, coupled with the high 
star formation rate in the galaxy, make it a good candidate TeV gamma-ray source. Observations were carried 
out in 2000 and 2001 with the CANGAROO-II 10 m imaging atmospheric Cerenkov telescope. TeV gamma-ray 
emission is detected at the ~ IIct level with a flux of (7.8 ± 2.5) x 10“'^^cm“^sec“^ at energies >0.5 TeV. The 
data indicate that the emission region is broader than the point spread function of our telescope. 
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1. Introduction 


NGC 253 is a nearby (~2.5 Mpc) (Vaucouleurs 1978) 
starburst galaxy, in which a high cosmic-ray density and 
non-thermal emissions are expected (Voelk et al. 19961 ). 
The large-scale radio continuum structure of NGC 253 
was studied by Hummel et al. (1984). Carilli et al. (1992) 
found a large, bright synchrotron-emitting halo extending 
to more than 10 kpc from the galaxy centre. Diffuse X-ray 
emission from the halo has also been detected (Cappi et 
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al. 199£; Strickland et al. 2002). The OSSE instrument on 
the Compton Gamma-Ray Observatory (CG RO) d etected 
sub-MeV gamma-rays (Bhattacharya et al. 1994| ) which 
were attributed to far-infrared photons scattered by the 
high energy, synchrotron emitting, electrons (Goldshmidt 
& Raphael! 1995| ). The EGRET detector, however, gave 
only upper limits on gamma-ray emission in GeV energies 
(Blom et al. 1999). 


We have observed NGC 253 with the CANGAROO-II 
telescope between October 3 and November 18, 2000, and 
September 20 to November 15, 2001. While the detection 
of diffuse TeV emission from our galaxy is very difficult 
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(e.g., Aharonian et al. 2002), observations of NGC 253 
provide the opportunity to study the distribution of cos¬ 
mic rays in galaxies like our own. 


2. Observation and Analysis 

Very high-energy gamma-rays (with energy greater than 
several hundred GeV) interact in the upper atmosphere 
and produce cascades of secondary photons and particles. 
Secondary electrons and positrons with energy greater 
than the (altitude dependent) Cerenkov threshold emit 
optical photons in directions approximately the same 
as that of the incident gam ma-ra y. The CANGAROO- 
II telescope (Kawachi et al. ^001 ) was constructed near 
Woomera, South Australia (136°47’E, 31°04’S) to detect 
the several nano-second duration flashes of Cerenkov light. 
The telescope consists of a 10 m reflector and a camera of 
552 pixels, each a 0.115° square photo-multiplier (PMT). 
From the detected Cerenkov light images, the directions 
and energies of the incident gamma-rays can be recon¬ 
structed. 

A total of ^75 hours of observations of NGC 253 
(“ON”) were made, with a similar amount of background 
(“OFF”) observations. Fach ON-source run tracked the 
center of NGC 253 through its highest elevations, with 
OFF-source observations offset in right ascension made 
in order to estimate the background. An edge-on galaxy, 
NGC 253 appears optically (based on red Digitized Sky 
Survey [DSS2] data; http://skyview.gsfc.nasa.gov/), 
as a disc ~0.4° by ~0.1° in size. This is larger than our 
angular resolution, but much smaller than the telescope’s 
field of view (FOV) of 1.84° x 1.84°. After selecting data 
taken at high elevation angles (> 70°) in good weather 
conditions, a total of 2959 min. ON-source and 2417 min. 
OFF-source data remained for further analysis. 

First, “cleaning” cuts were applied to the camera 
images, requiring pixel pulse-heights of greater than 
~3.3 photoelectrons, Cerenkov photon arrival times 
within 40 nsec, and clusters of at least four adjacent trig¬ 
gered pixels in each event. After these pre-selections, we 
carried out a shower image analysis using the standard set 
of image parameters, distance, length, width, and a (Hillas 


1985), combining length and width (after an initial dis¬ 


tance cut) t o assign likelihoods to each event (Fnomoto 
et al. 2002a , 2002b ). The likelihood for both a gamma- 
ray origin and cosmic-ray proton origin were calculated. 
The cut to reject background events was based on the ra¬ 
tio of these two likelihoods. After these cuts, the image 
orientation angle (a) was plotted. A gamma-ray signal 
appears as an excess at low a after the normalized OFF- 
source a distribution is subtracted from the ON-source 
distribution. As shown in Fig. pi, an excess of events with 
a < 30° is clearly observed for NGC 253. The number of 
gamma-ray-like events is 1652±149 (ll.ltr), correspond¬ 
ing to 0.56 gamma-rays per min. 

The energy scale (the conversion factor from ADC 
value to number of photoelectrons) was determined by 
adjusting the cosmic ray event rate and the relation be- 


1 80 
1 60 
140 
120 
100 
80 
60 
40 
20 
0 

0 20 40 60 800 20 40 60 80 
a(degrees) 

Fig. 1. Distributions of the image orientation angle (a): 
a) for NGC 253 and b) for the Crab. The points with er¬ 
ror bars were obtained by subtracting the normalized off- 
source data from the on-source data. The ratio of events 
in the higher a (> 30°) regions for on- and off-source 
data was used as the normalization factor, which agreed 
with the ratio of observation times to within 2%. The his¬ 
tograms were obtained from Monte Carlo simulations of 
gamma-rays from a point-source. 



tween the total ADC counts and the total number of hit 
pixels. This estimation also agreed with the calculated sin¬ 
gle trigger rate for pixel due to the Night Sky Background 
(NSB). We obtained the threshold energy for the gamma- 
ray detection to be ~ 500 GeV for a E~^-^ spectrum and 
^ 400 GeV for a E~^-^ spectrum, after pre-selection cuts 
and also the same values after image parameter cut. These 
values were obtained from the peak positions of the effec¬ 
tive area multiplied by power law energy dependences. 

The same analysis procedure was applied to Crab neb¬ 
ula data from observations in November and December 
2000 as a check. Approximately 10 hours of good data 
were obtained. The resulting a plot, shown in Fig. has 
an excess of 393±59 events (6.7cr). The solid histograms 
are the Monte Carlo predictions for the a distribution in 
the case of a point source. Note that the Crab observation 
was carried out at large zenith angles, around 56°. At such 
large zenith angles, the a distribution deteriorates due to 
the shrinkage of Cerenkov images. In order to compensate 
for this, we applied tighter cuts in the Crab analysis. As 
a result, the a distribution for Crab is slightly wider than 
that for NGC 253 (which is observed near the zenith). 
The experimental result, however, agrees with the Monte 
Carlo prediction in case of the Crab observation. For NGC 
253, the a distribution is broader than the point spread 
function (PSF) by a factor of two. Our Monte Carlo sim¬ 
ulations predicted that 73% of events from a point source 
would have a < 15°, whereas in fact 56% of events were 
contained within this range. 
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3. Results 

The spatial distribution of the gamma-rays from NGC 253 
was studied using the so-called signihcance map. The thick 
solid contours in Fig. ||are our results. The 65% confidence 
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Fig. 2. Prohle of the emission around NGC 253. The solid 
thick contours were obtained from our observations. This 
was made from the distribution of the detection signifi¬ 
cance determined at each location from the differences in 
the a plots (ON- minus OFF-source histogram) divided 
by the statistical errors. The solid thin contours are the 
DSS2 (red) data. The dotted contour was obtained with 
the CANGAROO-II telescope for the Crab. 


level contour, which is roughly one standard deviation of 
our angular resolution, is shown and compared with an op¬ 
tical image from DSS2, which is shown by the solid thin 
contours. The size of the TeV emission region is of the 
same order as the optical image, or larger. The dotted 
contour is the 65%-contour obtained by the Crab obser¬ 
vation. It is consistent with the Monte Carlo prediction of 
0.25°. That of the near-zenith observations was estimated 
to be 0.23°. In order to estimate the telescope pointing ac¬ 
curacy, we analyzed bright stars (optical magnitude from 
5 to 4) in other observations at various elevation angles, as 
the stars in the NGC 253 field of view are fainter than our 
detection limit. We conclude that the pointing uncertainty 
was within 0.1°, which is less than the pixel size. The 
center of gamma-ray emission from the Crab observation 
agreed within this accuracy. The gamma-ray acceptance 
of the telescope gradually decreases beyond 0.3° from the 
center of the field of view (dropping to half at 0.8°), and 
so we are unable to make a more detailed morphological 
study. It is clear, however, that the observed gamma-rays 
distribution from NGC 253 is inconsistent with emission 
from a point source. 

In order to derive the flux level, we estimated the ac¬ 
ceptances (effective areas) using Monte Carlo simulations 
(Enomoto et al. 2002aD . We compared our effective area 
with that of Whipple (Fig. 6 in Mohanty et al. 1998D , 
which are these after clustering and distance cut, i.e., be¬ 
fore image parameter cut. Our effective area agreed with 


Whippl e incl uding the energy dependences. According to 
Fegan ( 1996 ), the threshold of Whipple telescope was the 
same as ours. The sensitivity for the CANGAROO tele¬ 
scope was represented by that for Crab-type source, i.e., a 
point source with Crab intensity at zenith. We calculated 
its acceptance and assumed the same background level, 
with the cut of a < 15 degrees. It corresponds to ~8cr 
detection per hour, slightly higher than that of Whipple 
(Reynolds et al. 1993). A difference appeared after the 
shape parameter cut, i.e., likelihood analysis. Our accep¬ 
tance for the likelihood cut was as high as 86%. On the 
other hands, the “supercut” reduces it to 50% level (Fegan 
1996: Table II). The both background levels are the same. 
Note that our background level depends on NSB, artihcial 
light background and etc., because we used the experimen¬ 
tal OFF sample as the background sample for Probability 
Density Functions. 

The systematic uncertainty was calculated on a bin- 
by-bin basis by changing the cut values on the likelihood 
from 0.1 to 0.6 in steps of 0.1. The uncertainties are typ¬ 
ically ~30%. The energy scale uncertainties due to the 
mirror reflectivity, mirror segment distortions, and Mie 
scatterings were estimated to be 15% (point to point) and 
20% (overall). The energy resolution was estimated to be 
35% on an event by event basis. The energy resolution is 
bad because we could not assume source point event by 
event basis, i.e., we could not correct it using core dis¬ 
tance. The spill over effect, therefore, is considered to be 
large for such soft energy spectrum. In order to derive the 
spectrum of NGC 253 self-consistently, we adopted the it¬ 
eration method. At first we used an acceptance calculated 
by the simulation with a spectrum and we obtained 

the index of —3.7±0.3. We then iteratively used this htted 
value in the simulations to re-derive the spectrum. This 
process rapidly converged at an index of —3.75±0.27. The 
same iteration procedure was carried out using a function 
with cutoff. It also showed a good convergence. The dif¬ 
ferential fluxes obtained in 2000 and 2001 agreed within 
these errors. A check was made using Crab nebula data 
(with an energy threshold of ~2 TeV) and the derived 
flux agreed with previous measurement s with in the sta- 
tistic al err ors (~ 20%) (Tanimori et al. 199S ; Aharonian 
et al. 200C ). Also, the cosmic ray flux calculated from the 
background event rate was checked and found to match 
the expected rate within 10%. A more detailed descrip- 
tio n of th e analysis will be presented elsewhere (Itoh et 
al. 2002a ). The differential fluxes are plotted in the spec¬ 
tral energy distribution (SED) shown in Fig. ||, together 
with measurements from other energy bands (Carilli et al. 


1992; Hummel et al. 1984; Cappi et al. 1999; Strickland et 
al. 2002; Bhattacharya et al. 1994; Blom et al. 1999| ). 

A simple power-law fit to the TeV spectrum of 
NGC 253 deviates greatly from the EGRET measure¬ 
ments at GeV energies (Blom et al. 1999): a turn-over 
below the TeV region clearly exists. We fitted a power law 
with an exponential cut off, hxing the power-law index 
to —1.5, i.e., an inverse Compton-like spectrum (where 
the exponential term is exp(—-\/E/a), with a a free pa- 
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Fig. 3. Multi-band spectrum of NGC 253. The black 
squares were obtained by this experiment. When we con¬ 
strain th^flux to 3/4 of the EGRET upper limit (Blom 
et al. 1999 ) at 0.2 GeV using the function in the text, a 
flux proportional to £:-i-5g-\/^/(o.25±o.oi) obtained 
with x^/DOF= 1.8/5 (the dotted line). The X-ray data 
for CHANDRA and BeppoSAX soft component were cor¬ 
rected for photo-absorption in the Galaxy, but no correc¬ 
tion was applied to the BeppoSAX hard X-ray data. No 
correction was made for photo-absorption inside NGC 253 
itself. We note that the BeppoSAX, OSSE, and EGRET 
data were not able to spatially resolve NGC 253, in con¬ 
trast to the other data. 


rameter). The resulting improved fits are shown by the 
dotted curve in Fig. The minimum power-law index 
providing an acceptable fit (based on the value) was 
— 1.8 (2(7 limit). If the observed gamma-rays are produced 
by inverse Compton scattering of cosmic microwave back¬ 
ground (CMB) photons, then the fitted spectrum sug¬ 
gests that the maximum energy of the parent electrons 
is around several TeV. The integral flux is estimated to 
be (7.8 ± 2.5) x 10“^^cm“^s“^ at energies greater than 
0.52 TeV. 


4. Discussions 

The observations reported here are the first detection of 
TeV gamma-rays from a normal spiral galaxy (other than 
our own), and reveal the emission to be spatially extended 
and temporally steady. NGC 253 has been observed over 
a range of photon energies, as depicted in Fig. ||, however 
the SED is more difficult to interpret. The emission of 50- 
200 k eV photons observed by OSSE (Bhattacharya et al. 
1994 ) was interpreted as being inverse Compton scatter¬ 
ing of far-infrared photons from dust in the central region 
(Goldshmidt & Raphaeli 1995). A simple extrapolation to 
higher energies exceeds the EGRET 2a upper limit (Blom 
et al. 1999| ). Clearly, it is not possible to explain the OSSE 
and CANGAROO-II observations and the EGRET upper 
limits by invoking inverse Compton emission from a single 
population of electrons. The very large radio halo, extend¬ 
ing over ^10 kpc (Hummel et al. 1984; Carilli et al. 1992) 
and up to X-ray energies (Cappi et al. 1999), suggests 
the existence of a population of very high-energy cosmic 


rays which may be responsible for the inverse Compton 
production of TeV gamma-rays (Itoh et al. ^002b ), quite 
separate from sources concentrated near the centre of the 
galaxy which emit the majority of the photons observed 
by OSSE. 

NGC 253 is the first of a new class of object to be de¬ 
tected at TeV energies. The extended and steady nature of 
the emission makes it clear that the TeV gamma-rays are 
produced in a different environment than that of previ¬ 
ously reported extragalactic sources of the active galactic 
nucleus class. Studies of NGC 253 will, like those of the 
LMC at GeV energies (Sreekumar et al. 1992), provide 
the opportunity to learn more about the cosmic ray envi¬ 
ronment in galaxies like our own. 
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